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Abstract 

K^ ! In the context of anomaly-mediated supersymmetry breaking, it is natural 

ly-^ ! for vectorlike fields and singlets to have supersymmetry breaking masses 

^ ! of order 10 TeV, and therefore act as messengers of supersymmetry break- 

O '■ iiig. We show that this can give rise to phenomenologically viable spectra 



compatible with perturbative gauge coupling unification. The minimal mo- 
del interpolates continuously between pure anomaly mediation and gauge 
P-i! mediation with a messenger scale of order 10 TeV. It is also possible to 

Qh! have non-minimal models with more degenerate specta, with some squarks 

r-| ! lighter than sleptons. These models reduce to the MSSM at low energies 

and incorporate a natural solution of the ^ problem. The minimal model 
'k/< I has four continuous parameters and one discrete parameter (the number of 

H I messengers). The LEP Higgs mass bound can be satisfied in the minimal 

model by tuning parameters at the GUT scale to one part in 50. 



1 Introduction 

Anomaly-mediated supersymmetry breaking (AMSB) [1, 2] is an attractive mecha- 
nism for breaking supersymmetry (SUSY) without flavor problems. In this mech- 
anism, SUSY is broken by the VEV of a supergravity auxiliary field (-F^), whose 
couplings to matter are governed by scale covariance, and are hence naturally fiavor- 
blind. It defines a preferred renormalization group (RG) trajectory for all SUSY 
breaking couplings in terms of a single SUSY breaking scale (F^) ~ 10 TeV. Unfor- 
tunately, the slepton mass parameters are negative in the minimal supersymmetric 
standard model (MSSM). In this paper, we propose a solution to this problem based 
on an idea due to Nelson and Weiner [3], which built on early work by Pomarol and 
Rattazzi [4]. Nelson and Weiner considered a theory with extra vectorlike fields P 
and P and added a coupling of the form^ 

A£ = fd^e icPP + h.c. (1.1) 

This gives rise to a Dirac fermion mass c{F^) and scalar mass terms 

V = \c{F^)\\\P\' + |Pp) + (cl(F^)pPP + h.c.) . (1.2) 

The scalar mass-squared terms are positive for \c\ > 1. Assuming \c\ ~ 1, this is 
a supersymmetry breaking threshold at the scale (P,^), which gives SUSY breaking 
threshold corrections of order g'^ (F^) / 16tt'^ to SUSY breaking masses, taking them 
off the AMSB RG trajectory. As shown in Ref. [3], the leading threshold corrections 
to the scalar masses vanish, and the slepton mass-squared terms are therefore still 
negative at the scale (F^) . One can get positive slepton masses at the weak scale only 
by having a large number of messengers (5 or more 5 © 5's), which generates large 
gaugino masses at the messenger scale ~ 10 TeV, which in turn generates positive 
slepton masses from running between the messenger scale and the weak scale. How- 
ever, the resulting theories generally have charged slepton LSP, and the large number 
of messengers destroys perturbative unification. 

In this paper we consider a very simple extension of this model that has a more 
attractive phenomenology. The model consists of the MSSM plus a singlet S in 
addition to the vectorlike fields P, P. We include the most general interactions with 

^Couplings of this form with P and P replaced by the MSSM Higgs fields contribute to the 
Giudice-Masiero mechanism for generating the MSSM n term [5] . 



dimensionless coefficients. The additional terms in the Lagrangian are therefore 

A£ = fd'^e^ [icsS^ + cpPP) + h.c. 
"Ac 



+ d'd 



3! 



^S'^ + XpSPP 



(1.3) 



h.c. 



A superpotential couphng of the form SHuHd is assumed to be absent.^ For |c5| < 1 
the potential for S has a local maximum at S* = 0, so {S) ^ 0. This gives rise to a 
more general threshold with none of the problems of the minimal model. 



2 The Threshold 

In this section, we compute the SUSY breaking from the threshold. The scalar po- 
tential that arises from Eq. (1.3) is 



V 



csiFhS + IXsS' + XpPP ' + \cp{Fl) + \pS\\\P\^ + \P\ 



+ \{F^)\'acsS' + CpPP) +h.c 



(2.1) 



The potential is quadratic in P, P, so we look for a minimum with (P) = (P) = 0. 
In the appendix, we minimize the potential for real couplings and VEVs. We show 
that the global minimum preserves CP for 



cs <0 



and the we obtain 



{S) = -^[3cs + Jcsics 



El 
s 



2Xs 



-Cs + Jcs{cs 



This gives rise to a mass term for P, P that can be conveniently written as 

A£= fd'^e(f)MPP + h.c, 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



^For example, it may be forbidden by a discrete R symmetry S{9) h^ —S{i9), P{0) t-^ +P{i9), 
P{e) ^ +PiiO), Hu{iO) ^ +Huiie), Hdiie) ^ -Hdiie), u^ie) ^ -u'^iie), with all other fields 



where 

M = M[l + e^r{F^)l (2.6) 

In this parameterization r 7^ parameterizes the deviation from a supersymmetric 
threshold, i.e. r = gives a pure anomaly-mediated spectrum below the messenger 
scale. The model of Nelson and Weiner has r = —2. We then have 



M = cp(l+X)(F<^), (2.7) 

(2.8) 



2 + \x(cs + ^-Jcs{cs- 



X 



where 



^ _ MS) 



cs{F^) 2cp\ 



p^s 



(Scs + ^csics - 8)) . (2.9) 



This shows that all values of M and r are allowed, since 1 + X can be small and have 
either sign. (Note that this does not require any Yukawa couplings to be large.) In 
order to avoid a negative mass eigenvalue for the scalars P, P at the minimum, we 
require 

|(r + l)(F^)|<|M|. (2.10) 

We now evaluate the threshold contributions to the standard model fields due to 
the P fields. The general formulas can be obtained from the methods of Refs. [6, 7]. 
The soft SUSY breaking terms can be parameterized by higher superspace compo- 
nents of dimensionless couplings via 

^o = -44ln^, (2.11) 



l_d_ 
d_ 



^V2 = -T^^, (2.12) 



^^ = -27^^' (2-13) 



where all couplings are taken to be real superfields. In the present model, all SUSY 
breaking is contained in the conformal compensator and the P, P mass term, so we 
have 

W ^ 2^^^^ V ShO? " 9h^ ) • ^^-^^^ 



Note that this imphes the presence of mixed anomaly- and gauge-mediated terms for 
scalar masses, as first pointed out in Ref. [4]. In this way, we can obtain expressions 
for the soft masses at the scale M in the effective theory where P and P have been 
integrated out: 

<m = l(F,)^|-r^|^A' + 2r(r + 1)^^' - (r + 1)'^a}, (2.15) 

my^{M) = -{F^) [r/?; - (r + 1)/?,] , (2.16) 

A{M) = -j{F^) [rf3', - (r + l)/3,] . (2.17) 

Here primed (unprimed) quantities refer to the theory above (below) the scale M. 
The anomalous dimensions are defined by 

n dgi d\nZ 

A = -77] , 7 = -77] • (2.18) 

(7 m yU (7 m /i 

The expression for the scalar masses can be simplified in the case of fields with no 
Yukawa couplings to messengers, for which 7' = 7. We then have 

ml{M) = m^^^^sB + ^rir + 2){F^)'^Ap,, (2.19) 



where 



^OAMSB = -^(^0)'^A. (2.20) 

and Ap = jS' — 13. Similarly, we can write 

mi/2(M) = mi/2AMSB + -{F^)^(^g (2.21) 

A{M) = Aamsb - ^(i"0)A/3A. (2.22) 

These expressions explicitly display the fact that the soft masses reduce to the AMSB 
values in the limit r ^ 0. The scalar masses (but not gaugino masses and A terms) 
also reduce to their AMSB values for r -^ —2, as in the model of Nelson and Weiner. 
In the generalized model, all soft masses reduce to the gauge-mediated values in the 
limit r ^ 00 with r{F^) held fixed. For general r, the SUSY breaking spectrum in this 
model interpolates continuously between anomaly mediation and gauge mediation 
with a messenger scale of order 10 TeV (assuming all dimensionless couplings are 
order unity). 



As with the case of pure gauge- and anomaly-mediated SUSY breaking, Eqs. (2.15)- 
(2.17) are leading order results in a power series with subleading corrections sup- 
pressed by 0{{{F(f) /M'^Y) and 0{{r{F^) /M'^Y). In the present class of models, it is 
natural to have M ~ {F^),r{F^), where these effects may be important. They have 
been calculated for the case of pure gauge mediation, where they are known to be 
numerically small unless the SUSY breaking is tuned to be close to the instability 
limit F/M"^ — > 1 [8]. Because these corrections are UV finite, they do not depend on 
the regulator, and therefore depend on the conformal compensator only through the 
superfield mass of the messengers (see Eq. (2.5)). We can therefore use the results for 
gauge mediation with the replacement F/M"^ — > (r + 1){F^)/M'^. Since the stability 
limit is \{r + 1){F^)/M'^\ < 1 here as well, the corrections are small in the absence of 
fine tuning. 

3 The fx Problem 

In the context of AMSB, we cannot get a phenomenologically acceptable Higgsino 
mass by adding a fi term 



AC = fd^e fi(l)H^Hd + h.c. (3.1) 



since this gives rise to -B ~ (F^) ~ 10 TeV. One possibility is the NMSSM, where 
the VEV of a singlet gives the n term. However, it is nontrivial to get a negative 
mass-squared term for the singlet. Here we briefiy discuss another possibility within 
the MSSM that gives a more minimal model. 

We consider a mechanism originally proposed by Randall and Sundrum in Ref. [1]. 
We show that this mechanism can be made natural with appropriate broken symme- 
tries. We add a term to the Lagrangian of the form 



A£rs = Jd^9c{Y + Y^)^HuHd + h.c. (3.2) 



Here we have included factors of by canonically normalizing Hu^d but not the field 
Y. Expanding this out, we obtain the potential terms 



HuHd + h.c. 

(3.3) 



A£rs = [-cIF^WY + yt) + c(FjFy - h.c. 

'cFJiY + yt) + cFl] Jd^e H^Hd + h.c. 
We see that we can naturally get a vanishing B^ term at tree level if 

{Y + Y^) = Q (3.4) 



and all couplings and VEV's are real. This is natural by CP invariance, and we then 
obtain an effective n term 

/i = c{Fy). (3.5) 

The Bn term is generated from AMSB, giving rise to a model with only one additional 
parameter. 

It is crucial that the Y appears in the combination Y + Y'' . This is natural if the 
field Y is invariant under a shift symmetry 

Y^Y + iX (3.6) 

where A is a real constant. We must also forbid a term of the form 

A£= fd*ec''^H^Hd + h.c. (3.7) 

J (J) 

The discrete R symmetry 

Y{e) ^ -Y{te), Hu{e) ^ Hum, H^ie) ^ -H^ite) (3.8) 

forbids the unwanted term Eq. (3.7), and also has the feature that the lowest com- 
ponent of Y is odd, while Fy is even. ^ The VEV for Y that we need is therefore 
protected by this symmetry. In order to make the Yukawa couplings invariant, the 
standard model fields must also transform under the discrete symmetry, e.g. 

u%d) ^ -u'^iie), (3.9) 

with all other fields even. 

This shows that the term Eq. (3.2) with Y treated as a spurion provides a viable 
fi term in AMSB that is natural by symmetries. Effectively, it justifies the inclusion 
of a running n term into the AMSB RG trajectory. It does not explain why the n 
term is the same size as other SUSY breaking terms. We leave this for future work. 

4 Spectrum and Phenomenology 

We now discuss the SUSY breaking spectrum that results from this model. We assume 
that the messengers come in complete SU{5) multiplets, so that the gauge coupling 
unification in the MSSM is not an accident. The simplest possibility is then that 



■^This symmetry also forbids unwanted couplings between the singlet S and the Higgs fields if S 
transforms as S{9) i-^ —S{i6). 



the messengers consist of N copies of 5 © 5. For perturbative unification, we require 
A^ < 4. Under the standard model gauge group, these decompose into a doublet and 
a triplet, each of which can have different couplings cp and Xp (see Eq. (1.3)). These 
give rise to different values for r for the doublet and triplet messengers, and hence 
different SUSY breaking masses for colored and uncolored superpartners. We assume 
for simplicity that the A^ messengers have the same coupling {e.g. there can be an 
unbroken SU{N) symmetry in the messenger sector). This can be relaxed to obtain 
even more general spectra. 

For large r, the spectrum is close to that of gauge mediation. However, because 
SUSY breaking is driven by anomaly mediation, the gravitino mass is naturally of 
order (-F^), alleviating the gravitino problem. This may not be large enough for 
large r, but it is possible (and natural) to have masses for the gravitino and other 
gravitational moduli that are parametrically larger than (F^) with SUSY breaking 
dominated by anomaly mediation [9]. 

The simplest model is completely specified at high energies by M, F^, r2, r^, 
N, and fi. One parameter is eliminated by requiring that the Higgs VEV takes its 
experimentally determined value, so this model has four continuous and one discrete 
parameter."^ Of these, the dependence on the messenger scale is only logarithmic, 
since it just sets the scale for the RG running down to the weak scale. Explicit 
formulas for soft masses are presented in Appendix B. 

For illustration, the spectrum of superpartner masses at the messenger scale is 
shown in Fig. 1 as a function oi r = r2 = r^ for M = 50 TeV, for A^ = 1 and A^ = 4 
respectively. For r < we can obtain positive slepton mass-squared parameters, but 
the right-handed sleptons are lighter than the bino, giving rise to charged slepton 
LSP. We therefore focus our attention on r > 0. The spectra are still qualitatively 
similar to gauge- and anomaly mediation in the sense that colored superpartners are 
heavier than uncolored ones. For example, obtaining positive slepton mass-squared 
parameters requires r > 1, which then implies rriq > 5m^. 

Quite different possibilities exist if r2 7^ r^. In Fig. 2 we show an example spectrum 
with A^ = 1 and r^ = —1. We again require r > to avoid a slepton LSP. We see 
that the spectrum is more degenerate, and the SU{2)\y contribution to superpartner 
masses is comparable to SU{3)c- For r2 > 2, the superpartners charged under 
SU{2)w are the heaviest, followed by the gluino, then right-handed scalars and the 
Bino. Such spectra open up new regions of SUSY parameter space that may be 
interesting to explore. These spectra have a light stop, and therefore requires an 



^The top quark Yukawa coupling is fixed by demanding that the top quark mass has its measured 
value. 



Superpartner Masses at M=50 TeV with N=1 
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Superpartner Masses at M=50 TeV with N=4 
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Fig. 1. Spectrum of superpartner masses as a function of r = r2 = rs for M = 

50 TeV, and N = 1 (top) and N = 4 (bottom). For gaugino masses we plot \M\ and 
for scalar masses, we plot |m^|^/^ x sgn(m^). All masses are in units of ^^/(IGTr^). 



Superpartner Masses at M=50 TeV with N=1 and r3=-1 




Fig. 2. Spectrum of superpartner masses as a function of r2 for N = 1, M = 

50 TeV, and r^ = —1. For gaugino masses we plot \M\ and for scalar masses, we plot 
|m^|-^/^ X sgn(m^). All masses are in units of F(^/(167r^). 



additional contribution to the Higgs quartic. Possibilities include a "fat" Higgs [10] 
or large D terms from exotic gauge interactions [11]. 

We give some representative points in parameter space in Table 1, assuming r2 = 
r^ for simplicity. At the scale M we evaluate the soft-breaking parameters using 
Eqs. (2.15)-(2.17), and evolve them down using MSSM RG equations to the stop 
mass scale m^-. (Since we have small mixing in the stop sector, we simply use the 
common stop mass.) At the scale m^, we determine the /j, parameter by minimizing 
the one-loop effective potential. This includes the largest 2-loop corrections to the 
effective potential because we use a value of yt that includes 1-loop QCD corrections 
[12]. We then add by hand the 2-loop QCD threshold corrections to the higgs mass 
fn\o, although this is a small correction (< 2 GeV) for small stop mixing. 

The spectra given in Table 1 satisfy all experimental constraints. The most severe 
constraint is the LEP Higgs mass bound m/^o > 114.4 GeV. Because we do not have 
large mixing in the stop-sector, we require mi ~ 1 TeV to satisfy the Higgs mass 
bound, and the experimental constraints on the sleptons and LSP are easily satisfied. 
As we have large stop masses, these models are fine-tuned. 



Table 1. Sample MSSM spectra. All masses are in GeV. The main text gives the 
definition of fine-tuning. 





Point 1 


Point 2 


N 


1 


4 


r 


14.6 


6.45 


F, 


7.19 TeV 


6.34 TeV 


M 


201 TeV 


81 TeV 


/^ 


485 


425 


tan/5 


17.1 


17.7 


TTlhO 


115 


115 


^IL 


380 


330 


^R 


190 


150 


rriqL 


1220 


1170 


rriuR 


1170 


1130 


^dR 


1065 


1120 


m-ti 


1070 


1050 


m2 


1180 


1150 


nig 


880 


1280 


^x? 


80 


165 


Tuning 


170 


55 



We quantify the fine tuning by the sensitivity of the Higgs to varying parameters 
at the GUT scale. The Higgs mass is quadratically sensitive to the stop mass, but 
this is not a fundamental parameter in this model. The most sensitive fundamental 
parameter is 5'3(Mgut), so we define 



Fine tuning = 



93{M, 



GUTj 



dv 



(91nt; 



dgsiMcvT) dlngs{MGVT)' 



(4.1) 



Because the sensitivity is through the stop mass, the tuning increases quadratically 
with the stop mass, while the lightest Higgs mass increases only logarithmically. This 
means that the fine tuning increases exponentially as a function of the lightest Higgs 
mass. This phenomenon is intrinsic to the MSSM, not just the present model, and is 
illustrated in Fig. 3. Note that the fine-tuning is somewhat less for a large number 
of messengers, since QCD is non asymptotically free in this case, and therefore the 
sensitivity to 5'3(Mgut) is reduced. 
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Fine Tuning in g3(IVIGUT) 



500 

450 

400 

™ 350 

I 300 

^ 250 

~ 200 

150 

100 

50 



* g3 Tuning, N=3 
i g3 Tuning, N=4 




114 116 118 120 

m_h (GeV) 



122 



124 



Fig. 3. Fine-tuning in QsIMqut) as a function of lightest Higgs mass rriho for models 
with r > for A^ = 3 and 4. 



5 Conclusion 

We have constructed a well-motivated minimal model that naturally breaks SUSY in 
a flavor-blind way with a messenger scale near 10 TeV. The minimal model with one 
messenger has four continuous parameters and one discrete parameter, and can give 
rise to spectra that are very different from scenarios considered in the literature. These 
include "compact" spectra with colored superpartners close in mass to uncolored 
superpartners, a feature of the spectrum that may help with SUSY naturalness. 
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Appendix A: Minimization of the Potential 



We minimize Eq. (2.1) assuming that all couplings are real. It is useful to write the 
potential as 



v = xl 






\S\^+('-^] (^S^ + h.c. 
\ As J \2cs 



(A.l) 



and use units where cs{F^)/\s = 1. We see that the phase structure is completely 
determined by the dimensionless parameter 



cs 



Writing 

where s and 6 are real, we have 

V_ 

This is stationary in 6 for 



(S) = se^\ 



— = (1 + ^ cos 2e)s^ + s^ cos e + i/. 

Ac 



s = or sin^ = or cos^ = . 

4e 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



We consider these cases one at a time. 

The case sin^ = is equivalent to {S) = s = real. In that case, we find stationary 
points 



Consistency therefore requires ^ < |- It is easy to check that 

\/(s_)<\/(s+),\/(0) for e<0, 
V(0) <\/(s_),\/(s+) for 0<^<|. 



(A.6) 

(A.7) 
(A.8) 



It remains only to consider the third condition in Eq. (A.5). In this case, the 
stationary points are 



s = s± = ±2\ 



2^-1 



(A.9) 
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Reality of s and | cos^| < 1 are satisfied only if 

e>o. 



(A.IO) 



We have V^(s+) = V{s_), as we expect since CP is spontaneously broken. We can 
check that 



We conclude that 



where 



V{0) < V{s±) for ^ < 1, 
V{S±) < V{0) for ^ > 1. 



(s) = <; < e < 1 



costz-i- = =Fi 



e-1 



4^(2^-1)' 
Restoring the units, we obtain the formulas used in the main text. 



(A.ll) 
(A.12) 



(A. 13) 



(A. 14) 



Appendix B: Formulas for Soft Masses 

In this appendix, we give some explicit one-loop formulas for SUSY breaking 
The beta functions for the MSSM gauge couplings with N2 doublets and N3 
are 

b, 



A 



IGtt^ 



gl 



where 



63 = -3 + N3, 

62 = 1 + N2, 

h = ll + N2 + lN3. 



The one-loop anomalous dimensions are 

1 



7Q3 

7«3 



167r2 

1 
167r2 



fgi + ^gi + Igt - %* 
f ^3^ + fgl - %,^ 



masses, 
triplets 



(B.l) 

(B.2) 
(B.3) 
(B.4) 

(B.5) 
(B.6) 



13 



1 



7rf3 

IL 

7e 

iHu 

iHd 



167r2 

1 

167r2 

1 

1 
167r2 

1 
167r2 



'i9i + hi 



^gl + gl 



^al 



^gl + gl- Qvt 



Sgl + gf 



(B.7) 

(B.8) 

(B.9) 

(B.IO) 

(B.ll) 



For the quark fields of tlie first and second generation, tlie top Yukawa coupling 
contribution should be dropped. We do not include the other Yukawa couplings, since 
they are negligible for small tan/5. The beta function for the top Yukawa coupling is 

yt 



A 



yt 



167r2 



16^2 



13^2 



Qyt - fgt - 3^2 - fgl 



[B.12) 



These formulas can be used to compute the MSSM soft masses using Eqs. (2.19)- 
(2.22) in the main text. In the one-loop approximation, the contributions from the 
doublet and triplet messengers just add, and we obtain e.g. 

1 {F^? 



m 



Q,AMSB 2 mix 



2\2 



16^3-3^2-T^r + 2l/t(167r^A 



yt) 



(B.13) 



Am\ 



1 {F^f 



rs 



(rs + 2)N, {fgt + ^^4 



m. 



^ 2 (167r2)2 

+r,{r2 + 2)N2{3gt + lgt 
_ 1 {F^r 



««,AMSB 2 (167r2)2 

A 2 1 jF^r 

iXtti- = ^ — ^— — 

"« 2 (167r2)2 



16^3^ 



'-gl + Ayt{lQTi'(3, 



yt, 



[B.14) 
(B.15) 



f^3^3(r3 + 2)iV3 



+ f ^1 (1^3(^3 + 2)iV3 + r2(r2 + 2)N, 



m 



1 {F^? 



c(h,amsb 2 (167r2)2 
_ 1 {F^? 



-1^4 44 4 

16^3 - T^l 



(B.16) 
(B.17) 
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'^^ 2 (167r 



1 (i^^)^ 



2A2 



f(73V3(r3 + 2)iV3 
+ 1^1 (^3(r3 + 2)iV3 + r2(r2 + 2)iV2 



L,AMSB 2 ri67r 



2\2 



-3^2-ll^f 



[B.18) 
[B.19) 
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Am 



2 _ 1 W 



^ 2 (1671 



1 W 



2A2 



3gtr2{r2 + 2)N2 



+9t (1^3(^3 + 2)iV3 + r2(r2 + 2)iV2) 



^e«,AMSB - 2 n67r2)2 



-4%^ 



Am 



1 (i^^)^ 



m 



"« 2 (167r2)2 

_ 1 (F^)' 

H^AMSB - 2 (15^2)2 



4^1 (i^3(r3 + 2)iV3 + r2(r2 + 2)iV2 



-3gt-ngt + 6yt{167r^P. 



yt, 



[B.20) 
[B.21) 

(B.22) 

(B.23) 



m 



Hrf,AMSB 



2 167r 



1 {F4,f 



3(72V2(r2 + 2)Ar2 



+^t (1^3(^3 + 2)iV3 + r2(r2 + 2)N, 



Am 



2 (167r2)2 

1 {F^r 

"" 2 ri67r2)2 



-3^2-ll^t 
'igtr2{r2 + 2)Ar2 
+^f (1^3(^3 + 2)A^3 + r2(r2 + 2)Ar2 



[B.24) 

[B.25) 



[B.26) 



For the squarks of the first and second generation, we drop the top Yukawa couphng 
contribution. 



(B.27) 
(B.28) 
(B.29) 



The 


gaugino 


masses are 


given by 








mx^ 


167r2 V 


-ll + |r3iV 






mx2 


ic_2 V 


-l + r2iV2), 



TTlA, 



167r2 



(3 + r3iV3 



where the first term in the parenthesis is the AMSB contribution while the remaining 
terms are contributions from the doublet and triplet messengers. 
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